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1 Introduction 

While supermassive black holes probably gained most of their mass via accretion 
of gas, the galactic nuclei in which they are currently situated are dominated 
by stars. This article reviews recent theoretical work on the interaction between 
black holes and their stellar environment, and highlights ways in which the ob- 
served structure of galactic nuclei can be used to constrain the formation history 
of black holes. Nuclei may also contain dark matter, and the possibility of detect- 
ing supersymmetric particles via annihilation radiation from the Galactic center 
has generated some interest jHI ■ The evolution of the dark matter distribution 
in the presence of a black hole in a stellar nucleus is also discussed. 

2 Matter Distribution around Black Holes 

Figure^shows two natural models for the evolution of the stellar density around 
a black hole. If the star-star relaxation time is long, the form of p(r) will still 
reflect the way the nucleus formed. The simplest possible picture is growth of 
a black hole at a fixed location in space due to some spherically-symmetric 
accretion process. As the black hole grows it pulls in stars. If the growth time is 
long ( "adiabatic" ) compared with orbital periods, the change in p(r) is unique 
and straightforward to compute |51I64| . Figure QJi shows the result if the initial 
density is a power law, p oc r~'^° . The final profile is well described as two power 
laws joined at a radius rcusp where 

rcusp = arh, 0.19 < a < 0.22, 0.5 < 70 < 1.5 (1) 

and r/i is the radius containing a mass in stars equal to twice the final black hole 
mass. The final density within rcusp is 

pcxr"^, 2.29<7<2.5. (2) 

This result is insensitive to the geometry (flattened, triaxial, etc.) and degree of 
rotation of the nucleus. 

If the nucleus is older than roughly one relaxation time, a different sort of 
mass distribution is set up around the black hole, corresponding to a steady- 
state solution of the Fokker-Planck equation. This steady state has a nearly 
zero net flux of stars with respect to energy and the density profile is given 




Fig. 1. (a) "CoUisionless" density profiles, resulting from growtli of a black hole at a 
fixed location. Initial profiles were power laws, p on r , with 70 increasing upwards in 
steps of 0.25. The radial scale is normalized to Vh as defined in the initial (pre-black- 
hole) galaxy, (b) "CoUisional" density profiles, computed as steady-state solutions to 
the Fokker-Planck equation. The density inside ~ OAr^ satisfies p cx r~'^^^. 



uniquely as p oc r~^/^ near the black hole '57"3'. Figure^ shows the result, again 
assuming various power-law slopes for the initial p{r). The r~^/^ dependence is 
only reached at r < rcusp ~ O.lr/j. 

The condition for the appearance of a coUisional, 7 = 7/4 cusp is that the 
relaxation time, 



i3 



^G2mp(r) log yl' 

measured at ~ rh be less than the age of the nucleus |3I52| . Figure 01 shows 
Tr{rh) for a sample of elliptical galaxies drawn from Faber et al. Only a few 
galaxies, all low luminosity ellipticals, have < 10^° yr. Two examples are well 
known: the Milky Way bulge, which has « 2 x 10^ yr (assuming m = M©) , and 
M32, for which w 3 x 10^ yr. Both galaxies should have coUisionally relaxed 
nuclei. In the Milky Way, rh ~ 1.7 pc w 40" [T5| and rcusp ~ 4". The stellar 
density profile has been determined into ~ 0.1"; it is a broken power law with 
7 « 1.4 inside of ~ 10" and 7 « 2.0 outside ^ni- The break radius is reasonably 
close to the predicted value of rcusp but the interior slope is significantly flatter 
than 7/4. In the case of M32, r/, « 2.9 pc « 0.9" and rcusp ~ 0.1" f^. The 
stellar luminosity density has been determined into only ~ 0.15", not quite far 
enough to test for the existence of a coUisional cusp. The situation is worse for 
more distant galaxies that might harbor coUisionally relaxed nuclei, making the 
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Fig. 2. Dark matter density at t = (dotted lines) and 10^" yr (solid lines) after 
heating by stars in the Galactic bulge |38|. 

Milky Way nucleus the only system in which the theory can currently be tested. 
The fact that 7 < 7/4 in the Milky Way nucleus may be due to mass segregation 
which produces a flatter density profile for the less-massive stars |4I5()| . 

Particle dark matter, if present in galactic nuclei, would respond to the 
growth of the black hole according to the "coUisionless" model. Its density near 
the black hole could therefore be very high, resulting in a substantial rate of 
self-annihilations if the particles are supersymmetric neutralinos l^- The dark 
matter distribution would also evolve due to scattering off of stars, on a time 
scale close to T^, the star-star relaxation time |2 71381 Ifi] . In effect, the stars act 
as a heat source for the dark matter. FigureElillustrates this via time-dependent 
solutions to the Fokker-Planck equation. Given the appropriate normalization 
with respect to the dark matter density near the solar circle, this nuclear profile 
implies a potentially observable flux of annihilation products |38| . 

3 Feeding Rates 

In a spherical galaxy, stars or dark matter particles with angular momenta less 
than Jic K, ^/2GMtrt will pass within a distance rt of the center of the black 
hole. The "loss cone" is that set of orbits satisfying J < Jic{E). If we imagine 
that the black hole is embedded in a steady-state galaxy populated by stars with 
J < Jic, stars on loss cone orbits will be consumed in one orbital period or less, 
giving a transitory feeding rate of ~ Att^ f {E , 0) jf^{E) stars per unit energy per 
unit time. This is sometimes called the "full loss cone" consumption rate and 
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Fig. 3. Stellar tidal disruption rates in a sample of elliptical galaxies computed us- 
ing the Cohn-Kulsrud [7] loss cone boundary conditions |60l . Open circles are "core" 
galaxies and filled circles are "power-law" galaxies. Upper and lower dashed lines are 
disruption rates in power-law nuclei with p oc and r~^'^ respectively. 



it is an effective upper limit to how quickly stars can be supplied to the black 
hole. For instance, in a p ~ r^^ nucleus, the full loss cone feeding rate due to all 
stars with E > ${rh) is - a^rt/G^M, or ~ a^/Gc^ when n « = 2GM,/a^. 
This is a high enough rate to grow black holes of the correct size after 10^° yr, 
i.e. ~ IO-^Mq yr-i for a « 200 km s^i. 

Full loss cones have been postulated from time to time, perhaps first by 
Hills 122, later by P. Young and co-authors (for the "black-tide" model of quasar 
fuehng) |63I65| . and most recently by Zhao, Haehnelt & Rees IHH, who noted 
that the scaling is consistent with the M, — a relation. But while transient 
loss cone refilling during, say, galaxy mergers is imaginable, it is hard to see how 
a galaxy could arrange for the 10^ or more refilling events that would be needed 
to keep a loss cone continuously filled. 

Instead, the loss cone is probably usually empty at radii r < Vcrit , where rcrit 
is the radius at which the scattering time into the loss cone equals the orbital 
period. Below Vcrit, stars are lost the moment they enter the loss cone, while 
above rcrit, stars may scatter in and out of the loss cone in a single orbit. For 
M, Ri IO^Mq, rcrit ~ Tfi, and rcrit < Hi at larger M, 13 . The feeding rate 
due to stars inside rcrit is ~ M (rcrit) /Tr (rcrit) since in one relaxation time a 
star's angular momentum can change by of order itself. Outside rcrit, the loss 
cone is effectively full but the feeding rate drops off rapidly with radius since the 
capture sphere subtends such a small angle (the "pinhole" regime). One finds 
in a careful calculation that the flux into the black hole is strongly peaked at 
energies near '^(rcrit) |30I7| . 

Figure^shows tidal disruption rates computed in a sample of elliptical galax- 
ies, with rt set to the tidal radius for solar-mass stars (gOj. (For M, > 10^ Mq, 
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Tt < Tg and stars would be swallowed whole.) The numbers in this plot super- 
sede the estimates in which were based on the now-discredited Magorrian 
et al. PJI black hole masses. Disruption rates reach ~ 10~^ yr~^ in the faintest 
galaxies with verified black holes like M32. The dashed lines in that figure are 
the predicted rates in a singular isothermal sphere nucleus: 

iV « 2.5 X 10-3yr-i ( 1 -\ ( -^^] ' oc M.-"''^ (4) 

^ VlOOkms-V \W^MqJ ' ^ ' 

and in a nucleus with p oc r^'^/^: 

/ N 21/5 / , X -4/5 

VlOO km s~V \IQ^MqJ ^ ' 

where the latter expressions in each case use the Af, — a relation The tidal 
flaring rate is probably nearly independent of M, for M, < IO^Mq, and if dE 
galaxies contain black holes, they would dominate the total flaring rate |60j . 

A fraction 25% — 50% of a tidally disrupted star is expected to remain grav- 
itationally bound to the black hole 0. If this is the case, and if disruption 
rates have remained nearly constant over 10^° yr. Figure 13 suggests that stellar 
consumption could be a major contributor to the growth of black holes with 
M, < IOHIq. 

The feeding rates estimated above were based on a theory derived in the 
late 1970's to describe consumption rates of stars by massive black holes at the 
centers of globular clusters |3I3UI7| . Globular clusters are many relaxation times 
old, and the assumption of a coUisionally-relaxed steady state is built into the 
theory. For instance, the Cohn-Kulsrud 7 loss cone boundary conditions, used 
in Figure|21 are only valid in a coUisionally relaxed system. But relaxation times 
in galactic nuclei can be very long and this fact requires a re-thinking of loss 
cone theory. Figure 01 shows histograms of three time scales relevant to feeding 
rates. The bottom panel is the relaxation time (eq. OJl evaluated at the black 
hole's influence radius r^. exceeds 10^*^ yr in most galaxies. The middle and 
upper panels show two time scales more directly related to the establishment of 
a steady-state loss cone. The time scale for scattering of low angular momentum 
stars into the loss cone is 

Te{r) « 9{rfTr{r) (6) 

where 9{r) is the scattering angle as seen from a star at radius r; the square- 
root dependence of 6 on Tg reflects the fact that entry into the loss cone is 
a diffusive process. If the stellar phase space is initially well-populated by stars 
whose angular momenta are close to the critical value for capture Jic, encounters 
will set up a steady state near the loss cone in the scattering time corresponding 
to0 = Oi, » (rt/r)i/2, or 

Tic « -Tr. (7) 

r 

Figure 0] shows that Tic is shorter than 10^" yr in most galaxies. 
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Fig. 4. Time scales Tr, Ti^ and Tgap defined in the text, for the sample of elliptical 
galaxies in Fig. |3| Thick/thin histograms of Tgap are for mg/mi = (0.1,1), where 
m2/mi is the mass ratio of the binary black hole that preceded the current, single 
black hole. 



But suppose that the loss cone was initially emptied of all stars with peri- 
centers below some radius ro 3> r*; for instance, a binary black hole might have 
ejected such stars via the gravitational slingshot (§4). Once the two black holes 
had coalesced, the time to diffusively fill this angular momentum gap would be 

Tgap ~ —Tr » Tic- (8) 
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In the binary black hole model, most stars with pericenters 

Tp < Kau = (9) 

will have been ejected prior to the coalescence. Here ^ rt is the semimajor 
axis of the binary when it first becomes "hard," ^ is the reduced mass of the 
binary, and K Ki 2. The middle histogram in Figure^ shows Tgap for two values 
of the binary mass ratio, 7712/7711 = (0.1, 1). In every galaxy, 

Tic « Tgap < Tr. (10) 

Whereas Tic is generally shorter than lO^*' yr, Tgap exceeds 10^° yr in most 
galaxies and this is true for all "core" galaxies which show evidence of cusp 
destruction. 

A more careful calculation 61 , applying the time-dependent Fokker-Planck 
equation to the galaxies in Figure |3| shows that the time for the feeding 
rate to reach 1/2 of its steady-state value after loss cone depletion by a binary 
black hole is approximately 

with q — 7772/^1 the binary mass ratio. While this is a highly idealized model, 
it gives an indication of the time required for loss cone feeding to reach a steady 
state in galactic nuclei, and suggests that the tidal flaring rates in luminous 
"core" galaxies might be much lower than suggested by Figure 13 



4 Interaction of Binary Black Holes with Stars 

Binary supermassive black holes are inevitable consequences of galaxy mergers. 
Unless one of the black holes is very small, the two quickly form a bound pair 
at a separation a~ah = G/i/4(7^ « [/i/4(777i -I- 7772)]7'7i with ji the reduced mass. 
Subsequently the binary separation can decrease, but only if the binary is able to 
exchange angular momentum with stars or gas. (See |42j for a review of binary- 
gas interactions.) Stars that pass within a distance ~ 3a of the binary undergo a 
complex 3-body interaction followed by ejection of the star at a velocity of order 
Ii/mi2 times the binary orbital velocity (the "gravitational slingshot" ISD). Each 
ejected star carries away energy and angular momentum, causing the semi-major 
axis, eccentricity, orientation, and center-of-mass velocity of the binary to change 
and the local density of stars to drop |25I23I24I46I53I36I37| . 

In principle, A^-body studies can reveal both the long-term evolution of the 
binary as it interacts with stars, and the back-reaction of the binary on its stellar 
surroundings. But unless great care is taken, A^-body simulations are likely to 
give misleading results. The reason is that time scales for two-body scattering 
of stars into the binary's loss cone are shorter by factors of ^ A^/10^^ in the 
simulations than in real galaxies, and the long-term evolution of the binary is 
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Fig. 5. Physical regimes for the decay of massive black hole binaries |48|. 

likely to be dominated by spurious loss cone refilling, wandering of the binary, 
and other noise-driven effects. Values of N that are easily accessible to direct- 
summation iV-body simulation, N < 10^, are so small that stars tend to be 
scattered into the binary's loss cone at a faster rate than they are kicked out by 
the gravitational slingshot 0S| ■ This is analogous to the "full loss cone" around 
a single black hole and it guarantees that the binary will never run out of stars. 

iV-body studies are most reliable when characterizing the early stages of 
binary formation and decay. Due to algorithmic limitations, most such studies 
[ 91341331211613 51 have been based on galaxy models with unphysically large cores. 
The first simulations with realistically dense initial conditions ^39 47.4 3 showed 
that the stellar density around the binary drops very quickly after the binary 
first becomes hard, converting a p oc cusp to a shallower ~ cusp. 

But the long-term evolution of massive binaries in real galaxies, and the 
effect of binary evolution on nuclear structure, are not well understood. FigureEl 
summarizes the likely regimes. For small N, of the order easily accessible to 
A^-body simulations, coUisionality is so high that the binary's loss cone remains 
full. The binary's binding energy increases linearly with time, oa t, and the 
stellar density profile becomes gradually flatter. The same sort of evolution would 
occur in a triaxial nucleus where centrophilic orbits can continuously repopulate 
the loss cone lIBj. At the other extreme, in giant galaxies, relaxation times are 
very long and almost no stars would be scattered into the binary's loss cone. 
The decay stalls, leaving a "hole" in the stellar density [HZ]. A mechanism that 
might contribute to continued decay in this regime is re-ejection: stars ejected 
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Fig. 6. A'^-dependence of the decay rate of a massive binary in A'^-body simulations 
|58|. Each of the three runs started from the same initial state, an equal-mass binary 
with a binding energy of 10. For small A^, the long-term decay satisfies ~ t, the 
expected behavior in the full-loss-cone regime (Fig.l^. As A' increases, the lower degree 
of coUisionality allows the loss cone around the binary to remain partially unfilled and 
decay slows. 

with less than galactic escape velocity can continue to interact with the binary 
^8 . The most interesting, and least-understood, regime is the intermediate one, 
corresponding to galaxies in which two-body scattering contributes substantially 
to the binary's decay without completely re- filling its loss cone 66 . Reproducing 
this regime requires either A^-body simulations with N > 10^, or a hybrid scheme 
that combines a Fokker-Planck integrator with a central binary. FigureElshows a 
step in this direction: A^-body integrations for N up to 0.25 x 10® using a GRAPE- 
6 special-purpose computer, coupled with a chain regularization algorithm to 
follow close interactions of star and black hole particles [SHI ■ The A^-dependence 
of the decay rate is clearly visible. Similar results were obtained [313 using a 
lower-accuracy A^-body code. Increasing A'^ by an order of magnitude will soon 
be feasible and this will allow simulations in which the binary's evolution is 
dominated by slow diffusion into the loss cone, as in real galaxies. 

5 Brownian Motion 

Both single and binary black holes undergo a random walk in momentum space 
as their motion is perturbed by gravitational encounters with nearby stars. Sim- 
ple encounter theory predicts that a massive black hole reaches a state of ki- 
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netic energy equipartition with the stars, giving it a mean square velocity of 
{V'^) « 3(m/M,)(T^. This "gravitational Brownian motion" has been suggested 
as a mechanism for enhancing the supply of stars to a massive binary; indeed it 
has been asserted |54l()| that wandering of the binary's center of mass can itself 
guarantee a continued supply of stars. 

While Brownian motion probably does affect the decay rate of binaries in the 
A''-body simulations 0S| j it is doubtful that the effect is significant in real galax- 
ies. The Brownian velocity of single black holes is found in A'^-body simulations 
to be [221 

Im,{V') « ^ma' (12) 

where is the ID, mean square stellar velocity within a region r < 0.5r?i 
around the black hole and m is the stellar mass. Equation H12|l holds in galaxy 
models with a wide range of nuclear density slopes, black hole masses and particle 
numbers up to 10^, the largest values so far used in direct-summation A^-body 
simulations In the case of the Milky Way black hole, equation 1)12(1 implies 
Vrms ~ 0.17 km s~^ (assuming m = Mq), a little larger than the "equipartition" 
value. The reason is that the black hole responds to perturbations from stars 
whose velocities have themselves been increased by the presence of the black 
hole 

Brownian motion of a massive binary is larger than that of a single black 
hole, for two reasons [3^]. (1) Stars are ejected superelastically from the binary, 
imparting a greater momentum to the binary than they would to a single black 
hole. (2) The dynamical friction force acting on the binary is less than that acting 
on a single particle due to the randomization of the ejection angles. The first 
factor dominates and produces an extra fractional contribution to {V"^) that is 
approximately H / (32^/2^^ \n A) , with H 16 the dimensionless hardening rate 
of the binary |53| . Thus the Brownian motion of a binary is only slightly larger 
than that of a single black hole of the same mass. This prediction has been 
verified in TV-body simulations |47I35| . The rms displacement of a binary from 
its otherwise central location would be very small in a real galaxy, probably even 
less than the separation between the two black holes. 

The claim that wandering can guarantee a continued supply of stars to a 
binary is based on the following conceptual model j54l6) . As N increases, the 
wandering first drops as expected from equation H12(l . but at sufficiently large 
Mt/m the binary empties its loss cone, producing a low-density core in which the 
binary is easily displaced. The amplitude of the binary's center-of-mass motion 
increases, allowing it to interact with a larger number of stars. 

The motivation for this model was the complex A''-dependence observed in 
a set of binary decay simulations |54I6| : the decay rate first dropped with in- 
creasing N until iV « 2 X 10^, then remained constant when N was doubled. 
But the postulated iV-dependence of the wandering was never verified. Even if 
the interpretation in terms of wandering is correct, extrapolating the A^-body 
results from N < 10^ to N ^ 10^^ is problematic. Furthermore the galaxy mod- 
els were Plummer spheres which have a much lower degree of concentration than 
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real galaxies, and the wandering amplitudes inferred by these authors would be 
enormous if scaled to real galaxies, > 10^ pc, Extended simulations of binary 
hardening in realistic galaxy models with N > 10^ should soon settle this matter. 

6 Evidence for Cusp Destruction 

The structural parameters of elliptical galaxies display a smooth dependence 
on galaxy luminosity from dE galaxies up to the brighter ellipticals, My ~ 
—20 Above this luminosity, galaxies generally have cores, regions near the 
center where the density falls below the inward extrapolation of the outer profiles 
(although the density almost always continues to rise inward). It is natural to 
attribute the cores to stellar ejection by a binary black hole and to use the "mass 
deficit" as a quantitative test of the binary black hole model (Fig. [7|). In two 
studies |49l55j , pre-merger profiles for a sample of galaxies were constructed by 
extrapolating power laws, p cx r~^", inward of the core radius. Milosavljevic et 
al. jin] found mean mass deficits of ^ IM, for 70 = 1.5. More recently, Graham 
[T??] fit the outer luminosity data to Sersic profiles and defined the mass deficit 
with respect to this template; mass deficits were found to be slightly larger, 
about twice M,} These numbers are reasonable but it is difficult to say more 
given the difficulties cited above in interpreting A'^-body simulations of binary 
evolution. 

An interesting unsolved problem is the effect of multiple mergers on the 
central density profiles of galaxies containing black holes. A simple argument 
suggests that the net effect, e.g. the core size, should increase both with the final 
black hole mass and with the number of merger events. Consider the merger 
of two galaxies with steep central cusps and black holes of mass mi and 7/12; 
the result will be a galaxy with a shallower central profile and black hole of 
mass mi2 = rrii + m2 ■ Now imagine that that same galaxy is produced by the 
merger of two galaxies whose central density cusps had already been destroyed 
by previous merger events. The final merger would lower the central density 
still more, producing a larger core for the same mi2. Testing this prediction via 
A^-body simulations will be difficult for all the reasons discussed above, but it 
should be possible at least to test how the early evolution of binaries affects the 
central profiles of galaxies in multiple mergers. 

The idea that cores grow cumulatively in mergers was used by Volonteri et 
al. j59| to predict the evolution of luminous and dark matter cores in galaxies. 
In their most extreme ("core preservation") model, the binaries were assumed to 
eject as much stellar mass as needed to shrink to a separation where gravitational 
radiation losses induce coalescence. The result, after 10^^ yr, was cores with sizes 
of order 3r^. While intriguing, this model lacks a firm foundation in the A^-body 
simulations. 

^ Graham (2004) presents his mass deficits as significantly smaller than in the ear- 
lier study, but he bases his comparison exclusively on values of Mdef computed by 
Milosavljevic et al. using 70 = 2; the latter authors computed Mdef for a range of 
values 1.5 < 70 < 2.0. 




Fig. 7. Observed surface brightness profile of NGC 3348. The dashed line is a Sersic 
model fit to the large-radius data. The mass deficit is illustrated by the area designated 
"depleted zone" and the corresponding mass is roughly 3 x 10* Mq I19| . 



In small dense galaxies, a destroyed cusp would be expected to re-form via 
the collisional mechanism discussed above, if the relaxation time at the core 
radius is less than the elapsed time since the merger. Alternatively, steep cusps 
may be due to star formation that occurred after the most recent merger. 

There is one context in which formation of cores by black holes has been 
tested and shown to work. Globular clusters in the Large and Small Magellanic 
Clouds exhibit a range of ages and core sizes, and the mean core radius increases 
as the logarithm of the age JOI- The observed age dependence is reproduced 
very well by a model in which the first population of stars leave behind ^ IOMq 
remnants which sink to the center and displace the less massive stars before 
ejecting themselves via three-body interactions |4()j . This model may have some 
relevance to galactic nuclei if supermassive binaries stall long enough for third 
or fourth black holes to fall in. 



7 Black Hole Ejections 



If two black holes do manage to coalesce due to emission of gravitational radia- 
tion, linear momentum carried by the waves will impart a kick to the coalesced 
hole of amplitude ~ 10^^ c The recoil velocity depends in a complicated 
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Fig. 8. Effect on the nuclear density profile of black hole ejection. The initial galaxy 
model (black line) has a p ~ density cusp, (a) Impulsive removal of the black hole. 
Tick marks show the radius of the black hole's sphere of influence rn before ejection. A 
core forms with radius ~ 2rh- (b) Ejection at velocities less than escape velocity. The 
black hole has mass 0.003Mgai; the galaxy is initially spherical and the black hole's 
orbit remains nearly radial as it decays via dynamical friction. The arrow marks Th- 

way on the mass ratio and spins of the two black holes and on the spin ori- 
entations 22 ■ At present only the mass ratio dependence is well understood: 
the recoil velocity peaks at m2/mi ~ 0.38 and falls to zero at m2/mi = 
or 1. As discussed by Hughes, Favata & Holz in these proceedings, uncertainty 
about the strong-field behavior of the effect allows only plausible upper and 
lower limits to be placed on the magnitude of the kicks. However for moderately 
large spins and prograde capture, even the lower limits exceed 100 km s^^ for 
0.2 < TO2/TO1 < 0.6, and the upper limit estimates reach ~ 500 km s^^ for 
favorable mass ratios |44| . 

An important possible consequence of the kicks is ejection of black holes from 
galaxies, particularly in the early universe when potential wells were shallow 
^4 20 62 . Whether or not a black hole is completely ejected, displacement of 
the black hole transfers energy to the nucleus and lowers its density within a 
region of size ~ rh (Fig.lHjl. Impulsive removal at T4ick ^ Kjscapc produces a core 
of roughly constant density within a radius 2rh if the initial density profile is 
p cx r^^ . The net effect is greater if T4ick ^ 14scape since the black hole can return 
to the nucleus several times before settling down. For Vkick ^ 0.25T4scape the 
change in the density is negligible. Cores formed in this way have a similar size to 
those produced in the binary black hole model and complicate the interpretation 
of mass deficits. 
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